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ABSTRACT 

Aims. To investigate the properties of the nuclear molecular gas and to address the nature of the deeply buried source driving the IR 
emission of NGC 4418. 

Methods. IRAM 30m observations and basic non-LTE, single component radiative transport modelling of HNC, HCN, HCO + , CN, 
HC 3 N and H 2 CO 

Results. We find that NGC 4418 has a rich molecular chemistry - including unusually luminous HC3N 7=10-9, 16-15 and 25-24 
line emission - compared to the other high density tracers. We furthermore detect: ortho-H2CO 2-1, 3-2; CN 1-0, 2-1; HCO + , 1-0. 
3-2, HCN 3-2, HNC 1-0, 3-2 and tentatively OCS 12-1 1. The HCN, HCO + , H 2 CO and CN line emission can be fitted to densities of 
n = 5 X 10 4 - 10 5 cirT 3 and gas temperatures 7\=80-150 K. Both HNC and HC3N are, however, significantly more excited than the 
other species which requires higher gas densities - or radiative excitation through mid-IR pumping. The HCN line intensity is fainter 
than that of HCO + and HNC for the 3-2 transition, in contrast to previous findings for the 1-0 lines where the HCN emission is the 
most luminous. Assuming all line emission is emerging from the same gas, abundances of the observed species are estimated to be 
similar to each other within factors of 2-5. The most noteworthy being a high abundance of HC3N and a small-to-moderate abundance 
ratio between HCN and HCO + . 

Conclusions. We tentatively suggest that the observed molecular line emission is consistent with a young starburst, where the emission 
can be understood as emerging from dense, warm gas with an additional PDR component. We find that X-ray chemistry is not required 
to explain the observed mm line emission, including the HCN/HCO + 1-0 and 3-2 line ratios. The luminous HC3N line emission is an 
expected signature of dense, starforming gas. A deeply buried AGN can not be excluded, but its impact on the surrounding molecular 
medium is then suggested to be limited. However, detailed modelling of HC 3 N abundances in X-ray dominated regions (XDRs) 
should be carried out. The possibility of radiative excitation should also be further investigated. 

Key words, galaxies: evolution — galaxies: individual: NGC 4418 — galaxies: starburst — galaxies: active — radio lines: ISM — 
ISM: molecules 



' 1 . Introduction emission (Cagnoni et al. 1998) and a resolved radio continuum 

emission support the starburst interpretation. 

S ^?o! h , eir discove 7 h l th f Infrared Astronomical Satellite In contrasti s fert _ like infrared colors indicate that the dust 

03 ■ (IRAS), luminous infrared galaxies (LIRGs, L m > 10 1 L ) have fa NGC 441g ifj bd heated fe m obscmed AGN Imanishi et 

been studied at almost all wavelengths (e.g. Sanders & Mirabel a , (2004) find broad MR ^ ^ and a ^ HCN/HCO+ /=1 _ 

1996). They radiate most of their luminosity as dust thermal ^ ratio of >J _ whichi th ^ ±e nQtion of an 

emission in the infrared. However, the nature of the main power obscured AGN xhe lack of hard x m this context be 

toaSnVfbothri smfundel CtlC ^ 01 * ^ A ™ l ° the A ° N bei " g C ° mpt ° n ^ 

The molecular line ratio HCN/HCO + 1-0 has been suggested 

NGC 4418 is a nearly edge-on Sa-type galaxy with deep afj a diagnostic tool for identifying the effect of an AGN on its 

mid-infrared silicate absorption features suggesting that the in- sulTOUnding mo i eC ular ISM (e.g. Kohno et al. 2001, Imanishi 

ner region is enshrouded by large masses of warm (85 K) dust et al 2QQ4) The x _ rayfj from the AGN impact the interstel . 

(e.g. Spoon et al. 2001, Evans et al. 2003). Due to the high obscu- ^ medium (ISM) sunounding it> creating an XDR - an X-ray 

ration, it is difficult to determine the nature of the activity that is Dominated Region . since X -rays may penetrate deep into the 

driving the luminosity. The infrared luminosity to molecular gas surrounding IS M, the resulting XDR may potentially become 

mass ratio - L m / M(H 2 ) = 100 is high for a non-ULIRG galaxy large and affect molecular p r0 p er ties on large scales. Molecular 

(Sanders, Scoville, & Soifer 1991). Furthermore, NGC 4418 has Une ratios are therefore p 0te ntially useful as XDR indicators, 

a high infrared to radio continuum ratio q * 3 (q= the loga- According t0 models by e _ g _ Maloney et al. (1996), HCO + is 

nthmic ratio of far-infrared (FIR) to radio flux densities). This destroyed in XDRs _ while HCN may enjoy high abundances , 

FIR-excess may be caused by a young, somewhat synchrotron resulting in an elevated H CN/HCO + 1-0 line ratio. Recent XDR 

deficient, starburst (Roussel et al. 2003). The lack of hard X-ray modeK howevej - question the sugges ted underabundance of 

HCO + in XDRs (e.g. Meijerink and Spaans (2005), Meijerink 
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Table 1. Observational parameters and Gaussian line fits. 
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HC 3 N 10-9 


90.979 


27 


0.77 


0.8 ±0.08 


2110 


122 


6.4 


HC 3 N 16-15' 


145.561 


17 


0.70 


1.7 ±0.08 


2120 


130 


8.6 


HC 3 N 25-24 


227.419 


11 


0.53 


1.6 ±0.2 


2130 


140 


10.0 


HNCO 4-3 


87.925 


28 


0.77 


< 0.36 (3cr) 








H 2 CO 2-1 


140.839 


17 


0.70 


1.05 ±0.2 


2156 


145 


5.7 


H 2 CO 3-2 


225.698 


11 


0.53 


0.6 ± 0.3 


2140 


145 


3.6 


OCS 12-11 


145.947 


17 


0.70 


0.8 ±0.14 




174 


4.8 



a) Rest frequency, b) Beam width, c) Main beam efficiency d) Integrated intensity in T* K , uncorrected for beam efficiency and beam size. (7^ is 
transferred into the main beam brightness scale, T mb , through multiplying T* A with the beam efficiency rj mb .) For the 30m telescope, the relation 
between r mb and flux is 5 =4.95 Jy/K. e) Center velocity of the line, f) Gaussian fit line width, h) Line fit peak intensity. 

i) First spingroup, j) Including both first and second spingroup. k) Second spingroup, 1) HC 3 N 16-15 is contaminated by para-H 2 CO. Its contribu- 
tion is estimated to 20%. Fitted Gaussians and integrated intensities are only for the HC 3 N component (see text). 



similar - or HCO + to be somewhat overabundant. Thus, the 
interpretation of an elevated HCN/HCO + 1-0 line ratio is not 
straight forward, in particular when line opacities and excitation 
are taken into account. Furthermore, dense gas with low ioniza- 
tion levels may also have a relative underabundance of HCO + . 
Therefore, further observations and modelling is necessary to 
asses the impact of AGNs and starbursts on their surrounding 
ISM. 

We have carried out an initial molecular line search of 
NGC 4418 with the IRAM 30 telescope. The goal was to ob- 
serve multiple transitions of high density tracer molecules and to 
search for complex molecules to get a first handle on the molecu- 
lar ISM conditions of NGC 4418. High density tracer molecules 
require gas densities in excess of 10 4 cm -3 in order to be col- 
lisionally excited. This is two orders of magnitude greater than 
the 12 CO 1-0 line which is the most commonly used tracer of 
molecular gas in external galaxies. 

We find a surprisingly rich ISM chemistry in NGC 4418 - 
with very luminous HC3N line emission - and that the HCN 3-2 
line emission is fainter than both HCO + and HNC 3-2. In sec- 
tion 2, observations and results in terms of line intensities, ratios, 
and RADEX models are presented. In section 3, we discuss the 
interpretation of the luminous HC3N line and the other line in- 
tensities in the context of XDR and starburst scenarios. We ten- 
tatively suggest a starburst scenario and briefly present how this 
conclusion may be tested. 



2. Observations and results 

We have used the IRAM 30m telescope to observe 1mm, 2mm 
and 3mm lines of HNC, HCO + , HCN, CN, HC 3 N, Ortho-H 2 CO, 
CO and HNCO towards the center of the luminous galaxy 
NGC 4418. Observations were made in May and July 2006, and 
the system noise temperatures were typically 150 K and 500 K 
for the 3 mm and 1 mm lines, respectively. Pointing was checked 
regularly on nearby continuum sources and rms was found to 



vary between 1.5" and 2". Beamsizes and efficiencies are shown 
in Table 1 . 

First order baselines have been removed from all spectra. 
Weather conditions during both runs were close to normal sum- 
mer conditions with precipitable water vapour, pwv = 4.3- 
5.1 mm in May and pwv = 3.5—4. 1 mm in June. The 1mm ob- 
servations of HCN and HNC 3-2 were repeated in the July run 
to confirm the line ratios. We found that line intensities could 
be repeated within 15% suggesting that, in general, 1 mm con- 
ditions both in May and July were stable. We also repeated the 
H2CO 2-1 2mm line and we found that the line intensity was 
repeatable within 10%. 

2.1. Line intensities and ratios 

We have detected the following line emission: HNC 1-0, 3- 
2; HCO + 1-0, 3-2; HCN 3-2, CN 1-0, 2-1; HC 3 N 10-9, 16- 
15, 25-24; Ortho-H 2 CO 2-1 and 3-2; CO 1-0, 2-1 (see Fig. 
1, Table 1). We also tentatively detect the OCS 12-11 line at 
the edge of the HC 3 N 16-15 band. HNCO 4-3 was not detected 
down to an rms level of 2 mK. In Fig. 1 we have marked the posi- 
tion of the HC3N 7=10-9, V6=l-0 rotational-vibrational transi- 
tion, and the CH3OH 3-2 line. They are both 2<x detections if one 
considers the rms only. Further observations are required to con- 
firm their presence. Integrated line intensities and fitted gaussian 
line parameters are presented in Table 1. A para-H2CO 2-1 line 
is blended with the HC 3 N 16-15 line (para-H 2 CO is shifted +42 
MHz, or -82 kms 'from the HC3N line). A gaussian fit where 
we set the velocity of the line center to 2120 kms -1 and the line 
width to 130 km s -1 result in an estimate of the para-H2CO con- 
tribution to 20% of the total integrated line intensity. This im- 
plies an H2CO 2-1 ortho/para line ratio of >2 (see Table 1 for 
the Ortho 2-1 line intensity). In warm gas (7\;40K) the ortho- 
to-para ratio is close to the statistical ratio of 3 (e.g. Kahane et 
al. 1984). 

Line ratios, corrected for beam size and efficiency, are pre- 
sented in Table 2. When correcting for beam size we require an 
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Fig. 1. Spectra of detected species towards NGC 4418 (a: 12:26:54.63, 6: -00:52:39.6)(J2000). Note that the intensity scale is T* 
in mK, and that line intensities shown here are not corrected for beam size. Top three panels: Velocity resolution is 15 kms -1 . 
Left: HC3N 10-9 (position for 10-9 ve=l-0 rotational-vibrational transition indicated). Center: HC3N 16-15 - somewhat blended 
with para-H2CO (see text). Tentative line emission from OCS 12-11 towards the band edge (position of CH3OH 3-2 indicated in 
spectrum). Right: HC3N 25-24 (at the edge of the band the beginning of the CN 2-1 line (CN line intensity cannot be trusted). 
Second row of panels: Velocity resolution is 30 kms" 1 . Left: HNC 3-2 and 1-0, Center: HCO + 3-2 and 1-0, Right: HCN 3-2. 
Third row of panels: Left: CN 2-1 (note the two spin groups), Center: CN 1-0, Right: CO 1-0 and 2-1 . Last row of two panels: 
Left: Ortho H 2 CO 2-1, Right: Ortho H 2 CO 3-2. 



estimate on the source size 6 S . From Imanishi et al. (2004) one 
can deduce an upper limit to the HCN 1-0 source size of 2". We 
adopt this as # s for our calculations - assuming that all high den- 
sity tracer emission is emerging from the same region. Errors 



in line ratios are generally 20% including rms and calibration 
errors, but not including source size errors. 
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Table 2. Observed Line Ratios 2 



Species 


Transitions 


Line ratio 


ro 


?_1 /l_o 

Z. 1/1 U 




HCO+ 


3-2/1-0 


0.45 


HNC 


3-2/1-0 


0.80 


HCN 


3-2/1-0 


0.20 b 


CN 


2-1/1-0 


0.40 c 


H 2 CO 


3(1,2)-2(1,1)/2(1,2)-1(1,1) 


0.35 


HC 3 N 


16-15/10-9 


0.90 d 


HC 3 N 


25-24/10-9 


0.50 



a) If 9 S is the gaussian source size (which we have assumed to be 2" 
(see text)), and 9 b is the beam size for the line, the ratio of integrated 
intensities for transition X and ¥, corrected for beam and source-size, is 

>. , . fr mb (X)dv Jr*(x)rfv(fl2(x)+fl2), mbm 
assumed to scale according to: J . , = r , w , ; , — * — — Errors 

in line ratios are »20% including rms and calibration errors. 

b) : Peak to peak ratio is 0.1 - but integrated intensity ratio is 0.2. HCN 

1- data from Imanishi et al. (2004). 

c) : The integrated line ratio should include all spingroups for 1-0 and 

2- 1. We only have the brightest of the two main 1-0 spingroup in our 
beam since the other spingroup occurs +800 km s"'away from the main. 
Integrated intensity ratio is 0.4 - peak-to-peak is 0.3. d) Line ratio some- 
what uncertain due to para-H 2 CO contamination (see text). 



2.2. Radiative transport modelling 

We used the online mean escape probability (MEP) code 
RADEX (Schoier et al. 2005) for radiative transport modelling 
of the observed molecular line ratios. 

The online version of RADEX assumes spherical cloud ge- 
ometry. We assume that all line emission is emerging from the 
same gas. This is of course unlikely to be true, but will give a first 
handle on physical conditions and relative column densities. 

RADEX models the average condition for a single cloud 
type. Here, the emission has been modelled as emerging from 
one structure of linewidth 150 kms" 1 , dealing with the global 
line emission at once. A more realistic approach is to assume an 
ensemble of individual clouds with a filling factor. Furthermore, 
multiple cloud types, and/or gradients within clouds should also 
be considered in a more advanced model - including issues on 
self-gravitation, alternative excitation and dynamics. This will 
be presented in paper II, where we also will investigate the ra- 
diative excitation further. The aim of the simple approach here is 
to look at average relative abundances, and physical conditions 
indicated by the line ratios. 

2.2.1. HCN, HCO\ H 2 CO and CN 

The dust temperature of NGC 4418 is estimated to be high, 
rd=85 K - and emerging from a compact nuclear structure of 
diameter « 70 pc (Evans et al. 2003). The dense molecular gas 
appears also to be compact and at number densities exceeding 
10 4 cm" 3 , the gas temperature should begin to approach that 
of the dust. Thus assuming that the gas is also warm, we have 
varied the gas kinetic temperature from 7\=50 to 150 K in our 
RADEX models. We have varied the gas densities n within the 
range 10 4 - 10 6 cm" 3 . Lower densities should not result in sig- 
nificant line emission from any of the observed high density 
tracer species - unless affected by radiative excitation (see be- 
low). Higher densities are unlikely due to the moderate excita- 
tion. 



We searched for a solution that would both accomodate the 
greater HCN 1-0 luminosity than that of HCO + 1-0 (Imanishi et 
al. 2004) - and the reverse situation for the 3-2 line, and at the 
same time would fit the low excitation of HCN (and the higher 
excitation of HCO + ). Solutions within the constraints given 
above could be found for a density range of n = 5 x 10 4 - 10 s 
cm" 3 . Higher densities are not allowed since HCO + becomes 
thermalized. Lower densities give too low excitation even when 
the temperatures are allowed to exceed 150 K. The column den- 
sity ratio between HCN/HCO + varies between 3 and 5 for the 
whole range of solutions. For n = 5 x 10 4 cm" 3 the best fit is 
found for a temperature of 150 K, for n = 1 x 10 5 best fit is for 
T k =80 K. 

Ortho-H2CO is also moderately excited, and the 3-2/2-1 
line ratio agrees reasonably well with the above physical con- 
ditions, even if temperatures exceeding 80 K and 150 K, respec- 
tively, would give a better fit. The column density ratio between 
H2CO/HCO + ranges between 1 and 3 for all solutions. 

Online RADEX currently contains no datafiles for handling 
CN, but a rough comparison to CS can be made, and run 
RADEX for CS with the above model parameters which sug- 
gests that the CN line emission may also be fitted to the above 
conditions. CN appears to be on average equally abundant as 
HCO + . 

The 12 CO 2-1/1-0 line ratio is low (Table 2) and requires 
densities nt 10 3 cm" 3 and may not be fitted to the above model. 
This is either due to the 12 CO emission emerging from a low 
density nuclear component - or the notion of a 2" source size 
is not applicable for the lower transition 12 CO emission. OVRO 
observations by Dale et al. (2005) of the CO emission show a 
structure unresolved in their 5" beam. 

2.2.2. HNC and HC 3 N 

The HC3N and HNC line luminosities are surprisingly high 
— and significantly more excited than the other high density 
species. This points to either the possibility that HNC and HC3N 
are affected by radiative excitation - or that they are emerging 
from a denser and/or warmer cloud component than the other 
species. 

Collisional excitation: 

If we adhere to the same constraints on conditions and search 
in the same range as in section 2.2.1., we find that the HC3N 
7=10-9 and 16-15 lines can be fitted to a component of density 
n — 2 x 10 5 cm" 3 and temperature 7\=80 K. In this context, the 
HC3N column density is high, similar to the HCN abundance. 
The observed 25-24 line is, however, too bright to be fitted to 
this model. More than 90% of the 25-24 line would have to be 
emerging from another component with a density n > 10 6 cm -3 
to produce a luminous 25-24 line. Such a gas component would 
contaminate also the other high density gas tracer lines - where 
the low excitation does not leave much room for such a very 
dense gas phase. The filling factor of an = 10 6 cm" 3 gas com- 
ponent would also have to be low. 

A collisional exitation of HNC would require densities of 
nZ5 x 10 5 cm" 3 for the given temperature range - which in the 
contaxt of a one component model is not an allowed density 
range for the HCN, HCO + and H 2 CO line ratios. 

Radiative excitation: 
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Pumping of HNC by 2 1 .5 pm radiation via the degenerate bend- 
ing mode has been discussed by Aalto et al. (2007) as a possibil- 
ity to explain the large HNC/HCN 3-2 line ratio in NGC 4418, 
Arp 220 and Mrk 231. For HNC, a mid-infrared background 
with a brightness temperature of 50 K will start to successfully 
compete with collisions at gas densities n < 10 4 cirT 3 . Thus, 
significant contribution to the line emission may then originate 
from a gas phase of lower density than the critical density for the 
species. 

HC3N may also be pumped via a number of infrared bending 
transitions, e.g.(Wyrowski et al. 1999): 

- Bending mode v 5 =l-0 at 663 cirT 1 (956 K) with A„ ; = 2.2 
s _1 , 15 /an (HCC bend) 

- Bending mode v 6 =l-0 at 498 cnr 1 (718 K) with A ul = 0.15 
s- 1 ,20 y um (CCNbend) 

- Bending mode v 7 =l-0 at 223 cm" 1 (321 K) with A ul = 6 x 
10~ 4 s-\45 pm (CCC bend) 

In Fig. 1 we indicate the position of the 7=10-9 V6 rotational- 
vibrational line as a 2cr detection, but more observations are 
needed to show whether this feature is real, or just a baseline un- 
certainty. Rotational-vibrational transitions of HC3N have been 
observed in the Galaxy towards hot cores (e.g. Wyrowski et al. 
1999). 

Column density comparison 

If HC3N line emission is affected by radiative excitation, it 
may have a larger filling factor than the line emission from unaf- 
fected molecular lines, such as HCN. Column density estimates 
done under the assumption of collisional excitation from the 
same gas may therefore overestimate the relative HC3N abun- 
dance. 

For NGC 4418 this may largely be related to the mass frac- 
tion of lower density gas, in relation to the high density gas frac- 
tion. We expect the HC3N abundances to be lower in the low 
density gas phase since this phase should be more vulnerable to 
photodestruction. Furthermore, the mass fraction of low density 
gas in galactic nuclei is often found to be significantly lower than 
that of the high density gas mass fraction, even if the lower den- 
sity gas has a higher filling factor (e.g. Aalto et al. 1994, 1995). 
We therefore estimate the errors in the abundance estimate due 
to potential radiative excitation, as rather small: within factors 
of a few. Errors due to insufficient data, and model oversimplifi- 
cation may be significantly larger. Observations of more HC3N 
transitions will help establish to what degree radiative excitation 
may affect the line intensities of the rotational transitions. 

3. Discussion 

3.1. Interpreting the HC3N line emission 

The relative HC3N 10-9 line emission in NGC 4418 is one of the 
highest measured in a galaxy on these scales - comparable to 
that found for Arp 220 (Aalto et al. 2002; Martin-Pintade et al. 
in preparation). The implied HC3N column densities rival those 
of HCN suggesting unusually large abundances. 

How can we understand the luminous HC3N line emission? 
Below we list three scenarios and their expected impact on 
HC3N abundances: 

- Warm dense gas: Rodriguez-Franco et al. 1998, find 
large HC3N/CN abundance ratios in hot cores (warm, 
star-forming, dense gas) which they attribute to very large 
HC 3 N abundances (X[HC 3 N] « 10~ 8 ) caused largely by 
evaporation of icy grain mantles and low destuctive cosmic 



ray and UV fluxes. HC3N can also be formed via the 
neutral-neutral reaction: C2H2 + CN -> HC3N + H (e.g. 
Meier and Turner 2005). 

- PDRs: HC3N is easily destroyed by UV photons (destructon 
rate 5.5 x 10~ 9 s _1 ) and very rapidly via reactions with the 
ions C + and He + (Prasad & Huntress, 1980; Turner et al. 
1998). Thus, we do not expect luminous HC3N line emission 
from PDRs since it would quickly become destroyed by 
photodissociation or by reaction with C + . This is also 
observed by Rodriguez-Franco and collaborators towards 
PDRs in the Orion complex. 

- XDRs: X-rays can penetrate large columns of gas and 
are therefore not as good as UV photons at dissociating 
molecules. There are UV photons in an XDR, but mostly 
in the form of relatively weak, secondary far-UV photons. 
Their impact on HC3N destruction is likely less than in a 
PDR. An XDR is characterized (among other things) by high 
abundances of C + and in particular He + - that coexist with 
neutrals like CO and C. Both C + and He + are sources of de- 
struction of HC3N. The He + abundances in XDRs can be five 
times higher than those in PDRs according to to XDR/PDR 
models (Meijerink and Spaans 2005). In the absence of icy- 
mantle grain processing and of important ion-neutral foma- 
tion mechanisms, the neutral-neutral formation mechanism 
indicated above will be too slow to keep up with the destruc- 
tion proceeses. 

HC3 abundances are low in PDRs and we expect them also to be 
so in XDRs - although no direct obervations exist, and detailed 
modelling is required to determine abundances. From Galactic 
work, it is well known that the molecule thrives in hot cores and 
in cold dark clouds. It has been observed at high resolution in 
the nearby galaxy IC 342 by Meier and Turner (2005) where 
they find the 10-9 emission to be strongly correlated with 3 mm 
continuum emission interpreted as sites of on-going starforma- 
tion. 

3.2. The HCN/HCO + 1-0 line ratio and XDRs 

HCO + is suggested by e.g. Kohno et al. 2001 to be a tracer of 
the fraction of the dense gas which is involved in starformation. 
They compared the HCN/HCO + 1-0 line ratio in a selection of 
Seyfert and starburst galaxies and find that the relative HCO + 
1-0 luminosity is significantly higher in starbursts. According 
to models by Maloney et al. (1996) a deficiency of HCO + is 
expected near a hard X-ray source - in an X-ray dominated 
region (XDR) - resulting in an elevated observed HCN/HCO + 
line intensity ratio. Interestingly, elevated HCN/HCO + 1-0 line 
ratios are found in ULIRGs by Gracia-Carpio (2006) — and in 
general the line ratio appears to correlate with FIR luminosity. 
The authors attribute this to the presence of an AGN where 
the X-rays affect the chemistry to impact the HCN/HCO + 
abundance ratio. 

Abundances can not be directly deduced from the 
HCN/HCO + 1-0 line ratio - even if the emission is indeed 
emerging from the same gas. Excitation effects must be taken 
into account. Imanishi et al. (2004) report an HCN/HCO+ 1-0 
line ratio of 1.8 for NGC 4418, suggesting that this supports the 
notion of an XDR. We observe, however, the opposite case for 
the 3-2 transitions and we find that both the 1-0 and 3-2 line 
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ratios can be fitted to a rather moderate HCN/HCO + abundance 
ratio (see section 2.2.1). 

Furthermore, recent XDR modelling suggest that underabun- 
dant HCO + may not be an expected feature of an XDR. In 
contrast, for a large range of conditions, Meijerink and Spaans 
(2005) and Meijerink et al. (2006) find the HCO + abundances 
exceeding those of HCN in the XDR. In their models HCN is 
more abundant than HCO + at the edges of the XDR, but in these 
regions molecular abundance are generally low and do not con- 
tribute significantly to the total line luminosity. A discussion of 
this can also be found in Aalto et al. (2007). 

Relatively low HCO + abundances may be expected in young, 
synchrotron-deficient, starbursts which would also result in 
HCN/HCO + 1-0 line intensity ratios exceeding unity. An alter- 
native interpretation to an elevated HCN/HCO + 1-0 line ratio 
may therefore be that some Seyfert nuclei are surrounded by 
regions of very young starformation. If the starburst is young 
enough not to have produced a significant number of supernovae, 
then the relative HCO + 1-0 luminosity may be low compared to 
that of HCN 1-0. 

3.3. CNandHNC 

All XDR models predict high abundances of CN - a radical 
which is also abundant in PDRs (e.g. Rodriguez-Franco et al. 
1998). For both XDRs and PDRs the CN abundance is expected 
to be significantly higher (orders of magnitude) higher than that 
of HCN. For NGC 4418 CN is detected - but line emission ap- 
pears to suggest abundances less than HCN. 

Luminous HNC emission has been observed towards a sam- 
ple of warm galaxies (e.g. Aalto et al. 2002) and overluminous 
HNC 3-2 in three galaxies (Aalto et al. 2007). Low HNC abun- 
dances are observed towards hot, dense cores in the Galaxy (e.g. 
Schilke et al. 1992), while in a PDR, equal HCN and HNC 
abundances are expected (e.g. Aalto et al. 2002, Meijerink et al. 
2006). In dense XDRs {nl 10 5 cm 3 ) Meijerink et al. 2006 pre- 
dict HNC/HCN abundance ratios of two. Mid-IR pumping may 
boost the HNC luminosity relative to the HCN luminosity. For 
NGC 4418 the HNC emission may either emerge from a PDR 
component (and affected by radiative excitation) or it emerges 
from a deeply embedded, dense XDR. 

3.4. NGC 4418: Buried AGN or nascent starburst? 

We find a rich molecular chemistry in NGC 4418 and lumi- 
nous HC3N line emission which in our Galaxy is associated 
with warm, dense starforming cores. In this context, the ele- 
vated HCN/HCO+ 1-0 line ratio found by Imanishi et al. (2004) 
may instead indicate young starformation. As discussed in sec- 
tion 3.2, we therefore advice caution in using the HCN/HCO + 
1-0 line ratio as a diagnostic tool for XDR activity since its in- 
terpretation is ambiguous. It is clear that further XDR modelling 
is required, including a continued dicussion of the HCN/HCO + 
line emission ratio, as well as modelling of HC3N abundances in 
XDRs. 

For NGC 4418, we tentatively suggest that the line emission 
can be fitted to a combination of a warm, dense gas component 
not strongly affected by UV or cosmic rays (hot-core like chem- 
istry) supplying abundant HC3N and HCN — with an additional 
PDR component (providing CN, HCO + and HNC line emission). 
This scenario is consistent with the suggestion by Roussel et al. 
that the high ^-factor is due to a young starburst. They also pro- 
pose that NGC 4418 is probably in a somewhat later evolution- 



ary stage than other high-g objects since it does contain a syn- 
chrotron source, and contains both a nascent starburst and a more 
evolved component. 

3.4.1. Future observational tests 

Even though we suggest a starburst explanation for the observed 
line ratios, a deeply embedded AGN, surrounded by starforma- 
tion, cannot be excluded. Without high resolution information 
we cannot tell whether the HC3N line emission is emerging from 
further out than, for instance, the HCN or HCO + line emission. 
In general, observations at higher resolution will help resolve the 
dichotomy around XDRs vs. starformation. Below we list further 
observational tests of the molecular ISM: 

- Alternative, and more unambigious, XDR diagnostics are 
provided by highly excited (J > 10) rotational lines of 12 CO 
(Meijerink et al. 2006, Aalto et al. 2007). XDRs produce 
very warm 12 CO compared to, for example, PDRs. Other 
lines that can help distinguish between XDRs, PDRs and 
young starformation include HOC + , NO, CO + and CS (e.g. 
Meijerink et al. 2007) 

- Observe higher J lines of all detected molecular species. 
This allows for more detailed modelling of the physical con- 
ditions in the gas. For HC3N these observations - as well as 
a search for the rotational-vibrational lines - should reveal 
a possible pumping situation. This is also important for the 
HNC line emission. 

- High resolution observations of molecular lines to determine 
from where the line emission is emerging. The upcoming 
ALMA (Atacama Large Millimeter Array) observatory will 
provide unprecedented spatial resolution which will allow 
cloud-scale observations of many of the luminous galaxies. 
Dynamical studies of the nuclear gas will enable measure- 
ments of enclosed masses, and searches for molecular out- 
flows. 



4. Conclusions 

1 . We detect line emission from HNC, HCN, HCO+ /= 1-0 and 
3-2, CN 1-0 and 2-1, HC3N 10-9, 16-15, 25-24, Ortho- 
H 2 CO 3-2, 2-1 (and tentatively OCS 12 -11) towards the 
luminous, deeply obscured galaxy NGC 4418. All species 
have low excitation with the exception of HNC and HC3N 
which are more highly excited. 

2. The HC3N line emission is unusually luminous with the 10- 
9 line having similar intensity as CN 1-0. This suggests large 
abundances of HC3N, a molecule usually associated with 
young starformation. 

3. We find overluminous HNC 3-2 emission - where the HNC 
luminosity is a factor of 2.3 times that of HCN 3-2. This 
property led us to propose that HNC may be pumped by 20 
/im IR continuum (Aalto et al. 2007). The high excitation of 
HC3N is also consistent with radiative excitation. 

4. We find an HCN/HCO+ 3-2 ratio of about 0.5. Imanishi 
et al. (2004) report the opposite for the 1-0 ratio. The 
HCN, HCO + , H 2 CO line ratios can be fitted to a single 
cloud component with 7^=80 - 150 K and gas densities of 
5 x 10 4 - 1 x 10 5 cm 3 - and an HCN/HCO+ abundance ratio 
of *3. 

5. The HCN/HCO + 1-0 line ratio has previously been sug- 
gested to indicate the presence of an XDR associated with 
an AGN. We find, however, that the molecular line ratios are 
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also consistent with a young starburst - in particular the lu- 
minous HC3N line emission. We advice caution in using el- 
evated HCN/HCO + line ratios as XDR indicators since the 
interpretation is not unique. 
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